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Inspired by the models of A. Rebei and G. J. Parker^ and A. Rebei et. a/P, we study a physical 
model which describes the behaviour of magnetic moments in a ferromagnet. The magnetic moments 
are associated to 3d electrons which interact with conduction band electrons and with phonons. We 
study each interaction separately and then collect the results assuming that the electron-phonon 
interaction can be neglected. For the case of the spin-phonon interaction, we study the derivation 
of the equations of motion for the classical spin vector and find that the correct behaviour, as 
given by the Brown equation for the spin vector and the Bloch equation, using the results obtained 
by D. A. GaranirP for the average over fluctuations of the spin vector, can be obtained in the 
high temperature limit. At finite temperatures we show that the Markovian approximation for the 
fluctuations is not correct for time scales below some thermal correlation time tth- For the case 
of electrons we workout a perturbative expansion of the Feynman- Vernon functional. We find the 
expression for the random field correlation function. The composite model (as well as the individual 
models) is shown to satisfy a fluctuation-dissipation theorem for all temperature regimes if the 
behaviour of the coupling constants of the phonon-spin interaction remains unchanged with the 
temperature. The equations of motion are derived. 

PACS numbers: - 



INTRODUCTION 



The discovery of the GMR effect in 1988, for which 
Peter Grunberg and Albert Fert were awarded the No- 
bel Prize in Physics in 2007, motivated scientific research 
of the magnetization dynamics at the scale of nanome- 
ters and led to the birth of a new field of research now 
called spintronics. Spintronic devices have nanometer- 
scale sizes, can operate in high frequencies (~ 1GHz) and 
have a wide range of applications which go from the cre- 
ation of small dimension (< 1/im) microwave frequency 
generators to the improvement of magnetic storage de- 
vices. 

To successfully design these devices one needs to de- 
velop the theoretical comprehension of magnetization dy- 
namics at the appropriate scales. The complete under- 
standing of magnetization dynamics at the nano-scale 
can only be achieved by theorizing from first principles 
and that implies a quantum mechanical treatment. In 
particular, if one wants to describe a spin system far from 
equilibrium, one needs to use the methods of quantum 
open systems far from equilibrium, namely the Keldysh 4 
or the LindblacP formalisms. 

It has been showrP^ that the linear coupling interaction 
of a spin with a bosonic bath allows for the existence of 
white noise in the equation of motion which, under some 
particular conditions regarding the density of states of 
the bath, adopts the form of the Landau-Lifshitz-Gilbcrt- 
Brown equation. Also, it has been shown that if the 
spin vector satisfies a Landau-Lifshitz equation supple- 



mented with white noise, then the magnetic moment as 
the average over the fluctuations of the spin satisfies, in 
the limit of low temperatures, a Landau-Lifshitz equation 
and, in the limit of high temperatures, a Landau-Lifshitz- 
Bloch equatiorP. The collection of these results together 
with the known result of formal equivalence between 
Landau-Lifshitz and Landau-Lifshitz- Gilbert equations 
allows one to conclude that the interaction with phonons 
(or other bosonic bath satisfying certain conditions) can 
be responsible for the motion of magnetic moments as 
described by a Landau-Lifshitz-Bloch equation which, 
in fact, gives a good description of the physical situa- 
tion at high temperatures^. A quantum field theoretical 
treatment of the s-d interaction of conduction electrons 
and spina^ has shown that, in the semi-classical limit, 
the magnetization obeyed a generalized Landau-Lifshitz 
equation. 

The need of increasing the speed of storage of infor- 
mation in magnetic media and the limitations associated 
with the generation of magnetic field pulses by an elec- 
tric current require the research for ways of controlling 
the magnetization by other means than external mag- 
netic fields. In 1996, subpicosecond demagnetization in 
ferromagnetic Nickel was achieved using a 60 fs laser in 
the experiments of Beaurepaire et. aiP. Manipulating 
magnetization with ultrashort (of the order of the fem- 
tosecond) laser pulses is now a major research challenge 
because at such time scales it might be possible to reverse 
the magnetization faster than within half a precessional 
period 8 . Because of this, it is of fundamental importance 
to understand the time evolution of magnetic moments at 
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high temperatures and time scales which are approaching 
the femtosecond. Reviews on the state of art of ultrafast 
spin dynamics and prospects are given, for example, by 
A. Kyrilyuk et. aM and G. Zhang et. aP. 

The ultrashort laser pulses are expected to strongly 
couple lattice oscillations, that is, phonons, and/or con- 
duction electrons to the spins. This suggests that we may 
consider an effective microscopic theory for the system 
in which the fundamental interactions are spin-phonon, 
spin-electron and electron-phonon type. In this work we 
introduce and study a theory to model this physical sys- 
tem in which, for the sake of simplicity, we will neglect 
the electron-phonon interaction. 

This paper is divided in three major sections. In 
the first section we consider a generalized version of the 
model of A. Rebei and G. J. Parker-2, that modelled the 
interaction of a spin-j and a bosonic bath, now written 
to describe the interaction between a spin-j and a bath 
of phonons which are assumed to be spin-1. The second 
section considers a model of a ferromagnet, inspired by 
the work of A. Rebei et. aZP, assuming the interaction 
of a magnetic system associated to 3d electrons, repre- 
sented by a collection of spin-j vectors, with conduction 
4s electrons as the bath. In both sections we use a path 
integral representation for the coherent state matrix ele- 
ments of the reduced density matrix associated with the 
system. In the case of the phonons, since we consider a 
linear interaction, they can be exactly integrated out. Af- 
ter obtaining the effective action appearing in the phase 
of the path integral we use a stationary phase approx- 
imation and obtain an equation of motion for the clas- 
sical spin vector. We then introduce a random field in 
the equations of motion through a Hubbard-Stratonovich 
transformation in the path integral expression. The high 
temperature limit is discussed and the Landau-Lifshitz- 
Gilbert-Brown and Landau-Lifshitz-Bloch equations are 
recovered in this limit. The finite temperature case is 
considered. In the case of the conduction electrons, the 
same procedure of integrating out the bath is not possible 
due to the non linearity of the interaction. Nevertheless, 
we use a Hubbard-Stratonovich transformation and ob- 
tain an expansion for the effective action of the system 
which we truncate at the second order. The third section 
compiles the results of the first and the second sections 
in a single model assuming a quantum system composed 
of a spin-j vector field interacting with a bath of spin-| 
electrons and spin-1 phonons. 

Finally, we list our conclusions. 



I. SPIN-PHONON THEORY 

In this model the system considered is, for simplicity, a 
single spin-j which interacts, linearly, with a bath of spin- 
1 phonons which can have longitudinal and transverse 
polarizations. 

The total Hamiltonian is written in the form 



We now introduce the quantum operators associated to 
the spin-j particle S a which are Hermitian operators sat- 
isfying the angular-momentum commutation relations 



S a , Sp 



IE 



(2) 



where £ a ^ 7 is the Levi-Civita tensor density (repeated 
indexes are assumed to be summed from now on unless 
otherwise stated) . We also introduce the creation and de- 
struction operators associated with the phonon degrees 
of freedom a^(k), a^k) which are labelled by the polar- 
ization index \i = — 1, or 1 and by a momentum index 
k. They satisfy the Weyl algebra 



[a M (k),at(k')] = ^„<5 3 (k-k') 



(3) 



The spin operators and the phonon operators commute. 
The system Hamiltonian is considered to be 

Us = -S a H a = — S • H. (4) 

We write the phonon Hamiltonian in the form 

n R = J2^H( k M^ ( 5 ) 

where w M (k) are the energy eigenvalues (transverse os- 
cillations have the some frequency, that is, w_i(k) = 
Wl (k)). 

We write the most general form for an Hermitian 
Hamiltonian coupling linearly the spin, in the macrospin 
approximation^, and the phonon operators, 

& = " E (K^H(k)S a + H a ,(k)S a a^(k)) . (6) 

k 

The SU(2) invariance is explicit if we ensure that the 
operator 

H a = {K^H (k) + H OM (k)o^(k)) (7) 

k 

transforms as a vector under a 5J7(2)-transformation as- 
sociated with the spin-j and consequently Hi will behave 
as a scalar. We can then write the interaction Hamilto- 
nian as 



l~Li — — S a H a — — S • H 



(8) 



Following Rebei and Parker^, we now consider the re- 
duced density matrix associated with the system, that 



is, 



Ps(t) = Tr R {u{t;t )p(t )Ui{t;t )} 



(9) 



H = H S + Hi 



Hi 



(1) 



where U(t; t ) = Texp (^— i J t * y.(s)ds^j in which T is the 

time ordering operator. In many problems of magnetic 
systems one can consider that the bath relaxes much 
faster than the spin. Then the system and the bath, 
at time to, are decoupled. We make also the assumption 
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that, at time to, the bath is in equilibrium with the spin. 
With these considerations, we have 



p{to) = Ps(*o) ® Ph(*o), 



(10) 



p R (t ) = Z^ exp 1-pH 



where Z R = Tr |cxp P%Rj j and /3 denote, respec- 
tively, the partition function and the inverse of the tem- 
perature of the bath at the initial instant of time. 

We will work the path integral expression for the re- 
duced density matrix associated with the system using 
coherent state basis for the Hilbert spaces associated with 
our physical problem. For the phonons we use the so- 
called holomorphic representation in which the coherent 
states are defined by 



||a)=ncxp(at(kK(k)) |0), 



which satisfy the useful property 

a M (k)||a)=a M (k)|0) 



(11) 



(12) 



These states provide a decomposition of the identity (of 
the reservoir) of the form, 



(13) 



Ifl = / rf/u- (a;) 1 1 ct) (o; 1 1 , 



-Li 7T 



(k)a M (k) 



For the case of the spin-j we use the so called spin co- 
herent statef^K^, 



S} = (l + |C(S)| 2 r J exp(C(S)S_)bj>, (14) 



with 



C(S) = e^tan-, (15) 
S = (S a ) = j (sin 9 cos ip, sin 9 sin ip, cos 9) T , (16) 



recall that S± — Si ± iS 2 and that \jj) is the highest 
weight vector of the representation of the group SU(2) 
labelled by j. These states have the property, 



(SISJS) 



S 



(17) 



and they provide a decomposition of the (system) iden- 
tity operator in the form, 



i s = J d/x(S)|S) <S|, 
dp(S) = 2 -^6(S 2 -f)d 3 S. 

4-7T 



(18) 



The properties (12 1 and (17) allow us to make a corre- 



spondence between classical and quantum quantities. We 
will consider path integral representations which make 
use of coherent state matrix elements which can be com- 
puted trivially, by the prescription (a^(k), a^k), S a ) — ► 
(a* (k), a M (k), S a ). One can also define a holomorphic 
representation in the case of the spin coherent states 
and general rules for the insertion of operators can be 
derived^. 

The spin coherent state matrix elements of the reduced 
density matrix at time t, 

p s (S f ,S i ,t) = {S f \p s (t)\S i ) 
= (S/lTrfljw^ioH^Wt^^Jls,), (19) 

can be written in terms of those at time to by using the 
propagator J(Sf, S it t; S 2 , Si, t ), given by 

ps{Sf,Si,t) = J dp(Si)dp(S 2 )J(Sf 1 S i ,t;S2,S 1 ,to) 

X psr(Si,S a ,td). (20) 

The propagator has a path integral representation of the 
form 



J (S/, t; S2, Si, to) 



s, 



P M (Si) / P/z(S 2 )exp(iI[Si,S 2 ]).F[Si,S 2 ] 



(21) 



where 

^[Si, S 2 ] = 5wz[Si] — 5wz[S2] 
+S S [S 1 ]-S S [S 2 ], 

in which Swz is the Wess-Zumino action, 



5WZ 



[S] = j dt du n • (d u n x d t n) , 

Jtn JO 



(22) 
(23) 



where 

n(u,t) = (sin(u9)cos(p,sin(u9)sinip,cos(u9)) T (24) 



is a map which continuously deforms the constant curve 
given by n = (0, 0, 1) T (the north pole) into the curve 
n(t) = S(t)/j through a great circle, a geodesic, in S 2 . 
We see that the Wess-Zumino action, given by equation 
(23), gives the area enclosed by the path traced by the 



spin vector and the two great circles from the north pole 
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of the sphere to the endpoints of that path. We also have and J"[Si,S 2 ] is the Feynman- Vernon influence func- 

tional which can be expressed as 



Ss[S] = - / dt(S(t)\H s (t)\S(t)) = 

Jt n 



J 

Jt< 



to 

dt S(t) • H; (25) 

to 



J-[Si,S 2 ]= J d/i(oo)d/i(ai)d/i(a 2 ) fc(aS,t;ai,t |Si) x Z^k(a\,-il3;a2,0\0) x /^(a^a^io^), (26) 



where we have defined the kernel 

fe(a^,i/;ai,ij|S) 

X> 2 a exp 



x exp 



i I dt L 



(a* f a(tf) + a(t i ) t a i ^ 
(27) 



in which 

iL=i(^a-at^)-i[« Jl (at(t) ) a(t)) 
+H l (a t (i),aW,S(t))] 



(28) 



and 



(a(t)\\H R (t)\\a(t)) 



Hi = 



(a(t)\\l R \\a(t)) 
(S(t)|( a (0||^(t)||q(t))|S(t)) 
(a(t)\\i R \\a(t)) 



(29) 
(30) 



Note that, in the above formulas, we have adopted the 
vector notation a = (a M (k)) and a) = (a*(k)) The 
Feynman- Vernon functional can be computed exactly be- 
cause the integrals are all Gaussian. After some algebra, 
we obtain 



J"[Si, S 2 ] = exp <^ -i 



£ £ dhdhiS^ih) S 2 . a (t 1 ))a 3 G aP (t 1 -t 2 )a 3 ( |^^) 



(31) 



where, 



k,/j. 

(n(w„(k)) + 1) 6{t) + n(w„(k))0(-t) 
n(w M (k)) + 1 



«K( k )) W„ rw 



(32) 



I 

we know that the field S(i) is associated with the classical 
spin and D(t) is associated with quantum and thermal 
fluctuations. The Feynman- Vernon functional in terms 
of these fields is 



J"[S,D] = exp j-i jf* dt, J* dt 2 (VtScih) G a ^(t! - t 2 ) ( ^ ( ( * 2 2 ) ) } , (34) 

where, 

iG a0 (t) = £ ^(kJe-^W* ( tf x + ^("^ ) ^,(k). (35) 

I 

In the stationary phase approximation one considers the to be zero. The equations of motion which follow from 
variation of the action in the phase of the path integral 



in which n(u>) = (e' 3uj — 1) _1 is the Bose-Einstein dis- 
tribution function. If we introduce the fields associated 
with the Keldysh representation 



s a (t) \ _ ( l(S a>+ (t) + S a ,_(t)) 

D a (t) J { S a ,+ (t) - S a ,-(t) 



(33) 
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taking this procedure are 



D(t) = S(t) x — - + D(t) x 



<ys(t) 



£D(t)' 

^eff 



S(t) = S(t) x i 5 ^ + -D(t) x p 
w w <5D(i) 4 w <5S(t) 



where 



S cff [S,D] = / dtT>(t) -H.- Hog T[S,T>}. 

J to 



(36) 
(37) 

(38) 



We notice that there is no quadratic term involving 
two S(i) fields. By performing a Hubbard-Stratonovich 
transformation to the quadratic term in the action which 
couples two D(i) fields, as irtH we obtain a random field 
in the equations of motion, £(t). The real part of the 
correlation function of this random field is 

Re{<&(<i)&(«a)>} = Re{(U0)^i - *a)» 



-iw M (k)(ti— i 2 ) 



k./f 



x (l + 2n(w M (k)))fl- /SM (k)}.(39) 
The last expression can be re-written as 

/•GO 

Re{]T / du;5(a;-a; M (k))^(k)^(k)e-^( k )^-^ 



x coth 



^(k) \ 



} 



(40) 



Let us assume that H a/1 (k) = H a (u)^(\s.)). In this 
case, we can further write, interchanging the sum and 
the integral and using the property f{ujk)8{uik — u>) = 
f{uj)5{ujk — of the delta distribution, we find 



/>oo 

Re{/ duy2 S ("-^W)H*(u)Hf,(w)e 



-iuj(ti—t 2 ) 



or, equivalently, 

/•OO 

Re{ / P (u)H^(u)H p {u)e- iu ^- t ''> 
Jo 



■ <'"' 11 I — )} (42) 



where we have defined the density of states 



(43) 



In the high temperature limit, the last expression be- 
comes 

2k B T Re{ pH H '^ H ^ e-^-^du}. (44) 
Jo w 

If we consider a linear dispersion relation for the longi- 
tudinal and transverse phonons, 



u>o(k) = cik, w±i(k) = c t k, 



(45) 



and if we take the continuum limit for the bath, we arrive 
at the density of states 

k,/J 



d 3 k 



S(u) - w M (k)) 



V ( 1 



2^ 2 \cf c 3 t 



(46) 



where V is the volume of the reservoir and po = p(u>o) 
is the density of states evaluated at some frequency ujq 
taken as reference. This means that if we want to recover 
the Gaussian random field in the limit of high tempera- 
tures we must ensure that 



1 



Re{H*(u;)Hp(u)} <x - 
In particular, one must have, 



(47) 



(48) 



so that, 



Re{(£ Q (t)^(i'))} = S aP 2a G k B TS(t-t'). (49) 

The constant olq is a parameter which gives the intensity 
of the Gilbert damping term in the equation of motion. 
The last condition constraints the form of H a (w). If we 
define 



( ^PO \ 1/2 TT ( \ 

Z a = 2 W H a (uj). 

then the above condition reads, 

z* a zp + zj 3 z a = 25 a p, 

which means that 

z a = e^- 

and that the phases must satisfy 

7T 

<Pa = <PP + (271+1) -, 

for all a different from /3 and with n being an integer. 
The form of H a (uj) is, thus, constrained to be 



(50) 

(51) 
(52) 

(53) 



H a (u) 



KPo \ 



-1/2 



(54) 



if the familiar fluctuation-dissipation theorem is to be 
satisfied. 

Assuming the above dependence for H a (ui) for arbi- 
trary temperatures, we find that the noise correlation 
function is written in the form 



= acOaB / — wcoth — - cos cut 
Jo n V 2 / 

= a G S a p(p{t), (55) 



G 



where we have defined the generalized function (p 
¥>(*) = 



00 duj , //3w\ 

wcoth I — — I coswi, (56) 



which completely characterizes the noise correlation func- 
tion. The study of this generalized function or distribu- 
tion can be found, for instance, in the book of C. W. 
Gardiner and P. Zoller on Quantum NoisS^. The above 
integral can be expressed as 



<p(t) = 



coth 

00 7 
aw 

-00 2tt 



/3w 



COS Ljt 



we 



(57) 



The first integral is convergent and can be computed. 
The result is 



POO 


du) 


Jo 








2tt 





coth 



Pcj 
~2~ 



cosech 



cos wi 



(58) 



The other integral can be computed if we regularize it 
with an exponential cut off, 



du> 
2^ 



we 



dw 
1 X 2 -t 2 



k (A 2 + i 2 )' 
Formally we can, thus, write, 



^ (t) = ~2^ 2C ° SeCh2 (^)* 



(59) 



(60) 



The correlation exhibits a thermal correlation time aris- 
ing from the behaviour of (p in the limit \t\ — ¥ 00 which 
is 



<p(i) -> -2tt (k B T) exp (-2irk B T\t\ 
= -27r(fc B T) 2 expf-^ 

V T Th 



(61) 



where rxh = (2irk B T) . In the case of arbitrarily 
small temperatures, this correlation time becomes very 
large (in fact, the correlation function becomes inverse 
quadratic in \t\). 

It is interesting to notice that the long time behaviour 
of the correlations is negative. One must study the full 
expression for <p for finite cut off A to understand what 
is happening. We have 



dt Re{(£ a (t)&(0))} = 2a G k B T6 aP (62) 



and what really happens is that, in fact, the cut-off de- 
pendent positive term is usually larger than the asymp- 
totic negative term, except at zero temperature when the 
two effects cancel each other. 



This whole discussion can be simplified if one under- 



stands that the integral in the definition of ip, (56), can 
be written as 



ip{t) = a G k B T— coth [irk B Tt] , 



(63) 



but then one must take into account, as pointed by G. 
W. Ford and R. F. O' Connel 15 , that the correct formula 
for the derivative of the hyperbolic cotangent is, in fact, 



dcotha; , 2 _ r/ N 

; = — csch x + 2d(x). 

dx 



(64) 



The results are identical in both approaches. 

The existence of a thermal correlation time, which in 
standard units is r Th = h/k B T = (1.27 x 1CT 12 s)/T, im- 
plies that the approximation of a Markovian description 
for the random field is only valid for times scales which 
are longer than this one or in the limit of large tempera- 
tures. In the case of the problem of controlling magneti- 
zation using ultrashort laser pulses which end up exciting 
lattice oscillations it might not be a good approximation 
to consider the process as Markovian. Equation (62) is 



remarkable because it is a manifestation of a fluctuation- 
dissipation theorem and it is, in fact, valid in all temper- 
ature regimes. 

We also note that, with the assumptions we've done re- 
garding the constants i? Q (w), we will always recover the 
Landau-Lifshitz-Gilbcrt form of the equations of motion 
in the limit D — > 0. To sec this we just need to note that 



SS e S 



5T>(s) 



= I dt'G(s-t')S(t') 

D(s)=0 J to 



(65) 



where the matrix elements of G are 



-<u„(k)*0(_t). (66) 



xe 



If we use all the assumptions we made, the last expression 
becomes 



G afs (t) = 2a G 8 a pj t [8{t)]6{-t), 



which yields the result 

SS e s 



SD(s) 



a G S(s), 



(67) 



(68) 



D(s)=0 



so that the equation of motion for the spin becomes al- 
ways the Landau-Lifshitz-Gilbert equation in this limit. 
In the case of high temperatures we also recover the 
Landau-Lifshitz-Gilbert-Brown equation with the ran- 
dom field and due to the results of D. A. GaranirP we 
are able to recover the Landau-Lifshitz-Bloch equation 
for the average over fluctuations of the spin. To do that 
we neglect a term in the equation of motion which cou- 
ples two spin fields and a random field. This is valid in 
the limit of small fluctuations. 
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II. SPIN-ELECTRON THEORY 

In this model the system can be viewed as a set of spins 
in a lattice, each spin written in the form of a general 
spin-j representation of the SU(2) group, modelling 4d- 
type electrons in a magnetic medium. We consider the 
reservoir to be composed of conduction electrons. The 
interaction Hamiltonian is an s-d type interaction of the 
conduction electrons with the spin vector field. 
The model Hamiltonian we consider here is again of the 
form of ([T]), where 



Hi • Si - J <d&i ■ J H ^ °» ( 69 ) 



given in the momentum representation by 
^s = -E H (- k )'S(k) 



k 

2 ^ 

k 



J(k)S(-k) • S(k), 



(70) 



The other terms are given by 



H R 



k,Q 



e(k)cT(k)c Q (k) 



-A^s(-k)-S(k) 

k 

-£s(-k)-h(k), 



(71) 



where e(k) are the energies of the conduction electrons, 
h(k) denotes the magnetic field felt by the conduction 
electrons and s(k) denotes the composite operator 



S(k) = i Yl 4(k'-k)(<r) Q ^(k') 



(72) 



k',a,P 



which is the momentum space representation of the spin 
density operator of the conduction electrons. In the 
above equation, er = (<7i, 02, (T3) T , is the vector whose 
components are the Pauli matrices. The Latin indices, 
{i,j,k, ...}, refer to space indices and the Greek indices, 
{a,/3, ...}, refer to spin indices. Furthermore, the opera- 
tors c a (k), cT(k) are electron annihilation and creation 
operators, which satisfy the algebra, 

{c a (kr), c^(k 2 )} = S^S 3 ^ - k 2 ), (73) 
{c Q (ki),^(k 2 )} = {cl (k!)4(k a )} = 0. (74) 



We will now consider the path integral representation of 
the reduced density matrix associated with the system as 
we did in section [T] As so, we are going to use a basis of 
coherent-states for the system and the reservoir (in this 
case we have to use Grassmanian variables for fcrmions) , 
so that the states of the Hilbert space are written in the 
form, 



where, 



|S>®||7>, 



(75) 



|S)=niS(k)>=n^(S(k))|0), (76) 

k k 

Il7> = LI l7«(k)> = Il ex P (4(k)7«(k))) |0) . (77) 



a,k 



a , k 



These states |S) are to be defined in such a way that 
they satisfy (S| S(k) |S) = S(k) (S| l s |S). Let us now 
consider the vector 

I s ') = + IC(S0| 2 r J eMt(Si)S-,i) Ho) , (78) 



where \ipo) is the tensor product of highest weight 
vectors of the spin-j representation and C( s i) = 
tan(6*i/2)e' Vi denotes the stereographic projection of 
S{ = j(s'm8 i cosipi,sm9 i smip i ,cos6i) T through the 
North Pole. Clearly, this state satisfies 



(s'|s i |s'> = s i (s , |i s |s , > = s i 



(79) 



we can now replace, in the last equation, Si and Si by 
their Fourier representations, yielding 

(S'|^e 4k ' x *S(k) |S') =^V k ' x *S(k), (80) 



5> 4k ' x *(S'|(s(k)-S(k)) |S')=0, 



(81) 



Since this relation is valid for all Xj we identify |S') with 
|S) because, indeed, one has 



<S'| S(k) |S') ^ S(k) <S'| i s |S') 



(82) 



Assuming that at time to equation (11) holds, we easily 



find, in analogy with section [I] th at the expression for 
the propagator is the same as |21[) but now in /[Si, S 2 ] 
one must sum over all spin-j degrees of freedom and the 
Feynman- Vernon functional is now given by 



^[Si,S 2 ] = / ^(70)^(71)^(72) fc(-7cU7i,*o|Si) x Z i , 1 fc( 7 I,-z/3; 72 ,0|0) x fc*( 7( U72,*o|S 2 ), (83) 



where d/x( 7 ) = JJ a k d 2 7a (k) exp (— 7 *(k) 7a (k)) and the in ( p7| (taking special care because Grassmannian vari- 
kernel fe( 7 t, t/; Ji, ij|S) is defined in the same way as 



ables anti-commute) . The last expression assumes, as in 
|IJ that the density matrix decouples at an initial time 
to and that the reduced density matrix of the bath, at 
that time, takes the form of a Boltzmann factor with 
Hamiltonian Hr and inverse temperature /3. The minus 
sign in the first kernel in ( 83 1 is due to the anti-periodic 



conditions on the trace formula using fermionic coherent 



states. Unlike the case of the phonons we cannot compute 
the Gaussian integrals exactly. We would then like to do 
some expansion depending on the parameter A. In order 
to achieve that we will perform a Hubbard-Stratonovich 
transformation and then expand a determinant resulting 
from the functional integrals. To do this, one can define 
an auxiliary bilinear form 



1 



( G ^) aP (* - k - k '; s ) = - 2 5{f ~ ^ [°" Q/3 ' h(k ~ k ' } + X(Tap ' s( *' k " k ' } ] ■ 

With this definition, it is clear that we can do a Hubbard-Stratonovich transformation of the form 

exp(-i [ f dsds' 7;(s,k)(G- 1 )^( S - S ',k-k';S) 7 , 3 ( S ',k')) 

J to J t kk',O0 



n^^exp^/' fdsds' J2 Cfek)(G)^( S - S ' ) k-k';S)^( S / ,k') + 
aet [&j j t0 j t0 kk%Q ^ 

/ / dsJ2C(s, kfrafa k) + j* a (s, k)Us, k)). 

Jtn Jtn 



+1 



(84) 

(85) 
(86) 
(87) 



k. • 



Doing this, we achieve a linear coupling between the 7's integrals in the 7's. Replacing this transformation in the 
and the £'s, which allows us to compute the Gaussian expression for the influence functional and computing the 

Gaussian integrals associated with the 7's, we arrive at, 



where, 



in which, 



det(-tG)[Si] det(*G)[S 3 ] 



x exp 



[i f fdsds' «?>> k ) a a (s,k)) (A(s-s',k-k') a0 ) ( ) I 



, Vik , 1 (an) a/3 (t,k) + (G)^(t,k;S 1 ) (Gi2) aP (t,k) 

; ~ ' (GnU (*. k ) (*. k ) - ( G W (*> k ; 



(Gu) a0 (t,k) (Gu) a p(t,k) 

(Q2i) aP (t,k) (g 2 2) a p(t,k) 

~ ld WW y l-/(e(k)) (l-/(e(k)))0(-t)-/(c(k))0(t) 

with f{x) = (e x + 1) _1 being the Fermi-Dirac distribution. 



(89) 

(90) 

(91) 
(92) 



The functional integral is readily evaluated to be, 



where, 



det{ 



Gn + GiSi) 
G21 



f G- 1 (S 1 ) 
I -G- 1 (S 2 ) 



012 
022 — G(S 2 ) 



} = det(l + 6), (93) 



SnG-HSi) -0i 2 G^ 1 (S 2 ) 
1 021G-HS1) -0 22 G- 1 (S 2 ) 



(94) 
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The determinant of equation ( 93 1 should be understood, 



of course, in the functional sense. We can write, 



det(l + 6) = exp{Tr [log(l + 6)]} 



here the trace is also understood in the functional sense. 
In the way it is written, this produces an expansion in 
powers of the matrix elements of & and consequently an 
expansion in the parameter A. 

The first term of the expansion is easily found to be 
zero and if we keep only terms of second order in the 
matrix elements of & we obtain, in the Keldysh repre- 



(95) sentation of the fields as defined by equation (331 



log [det (1 + 6)] = -2 Tr(32 + °( &3 ) 



X H I I d«idt2(«5ff)(*i-t2,k)e(ti-t 2 )h(-k)-D(t 3 ,k) 
X 2 J2 [ [ dt 1 dt 2 (Sg)(t 1 -t 2l k)0(t 1 -t 2 )S(t ll -k)--D(t 2 ,k) 
TE / / dtidt 2 (Vg){t 1 -t 2 ,k)n(t 1 ,-k)-T>(t 2 ,k)+0(& 3 ), 

t. J to J to 



(96) 



where 



and 



k' 



x(l-/(e(k')))/(e(k'-k)), (97) 



(V4>) (k) 
(S0) (k) 



1 



(0(k) + 0(-k)), 
(#k)-#-k)), 



(98) 



are the symmetrizer and the anti-symmetrizer operators 
associated with the momentum variable, respectively. 
This makes it clear that the part leading to dissipation 
will be the term coupling two D fields since it will intro- 
duce an imaginary term in the action. This is easy to see 



because the Fourier transform of an even (odd) function 
is real (imaginary), if the function is real. 

One can, analogously to what we did in section[TJ intro- 
duce a random field using a Hubbard-Stratonovich trans- 
formation. This time the correlation function (in position 
space) reads 

(£ Q (i,x)^(0)) = 6 aP ^-[(T- 1 oV)g](t,x) 

= <5^yRe{[J^- 1 9 ](t,x)}, (99) 

where T~ x denotes the inverse Fourier transform as- 
sociated with the variable k. We now prove that a 
fluctuation-dissipation theorem holds for this effective 
theory. In order to do that, we integrate the correlation 
function of equation ( 99 ) , 



ds [ d 3 x(^ a {s, X )^(0)) = S aP ^-Re{ [ ds f d 3 a:[J- 1 5 ](s,x)}. (100) 
Now we rewrite the integral appearing in the above expression by replacing the expression for g(t,h) of (97), 

f ds [ d*x[F- l g]{s,*) = f ds f ^^[k-x-^kV^k'-k)).] 

Jto J Jto J 

x(l-/(e(k')))/(e(k'-k)) (101) 



Identifying the Dirac delta's appearing in the above ex- 
pression, we find 



£<5 3 (k)<5( £ (k')-6(k'-k)) 



kk' 



x(l-/(e(k')))/(e(k'-k)), 



(102) 



or, if we define the excitation energy w(k,k) = e(k' — 
k)-e(k'), 

5> 3 (k)5( W (k',k)) 

kk' 

x(l-/(e(k')))/(e(k')+w(k',k)), (103) 
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The above formula has to be treated carefully. First we 
recall the useful identity, 

/(e)(l - f(e - w)) = »(w)(/(e - w) - /(e)), (104) 

where n(w) is the Bose-Einstein distribution, yielding 

£«J 3 (k)«$(a;(k',k)) 

kk' 

Xn(u,(k', k)) (/(e(k')) - /(e(k') + c(k', k))) , (105) 

Because of the delta associated with the excitation en- 
ergy we only need the integrand evaluated at w(k, k') = 
0. Clearly, the expression is undetermined because the 
Bose-Einstein distribution diverges and the difference of 
Fermi-Dirac distributions goes to zero. To proceed we 
will expand the integrand for small w(k, k'). Observing 
that /'(e) = — /(e) (1 — /(e)) we can Taylor expand the 
part containing Fermi-Dirac distributions so that we ob- 
tain 

n( W (k',k))(/(e(k'))-/(e(k')+^k',k))) 

= M^^ (k '' k)/(£(k ' ))(l " /(e(k ' )))] 
+0(u(k',k)) 
1 



The final result is 



(3 



/(e(k'))(l-/(e(k')))+G( W (k',k)). (106) 



This means that our sum becomes, simply, 



-y 



5 3 (k)^( w (k',k)) 



x/(e(k'))(l-/(e(k'))). 



(107) 



If we sum in k' the only contribution of the sum will be 
that in which w(k', k) = 0, or e(k' — k) = e(k). Assuming 
that e(k) is a power law in |k| this implies that the only 
term surviving is the one with k' = 0. This reasoning 
gives, 



i/(e(0))(l-/(e(0)))]T£ 3 (k) 

P k 

I/( e(0 ))(l-/(e(Q))) = ^. 



(108) 



ds / d 3 x(£a(s,x)^(0)) 



ds J d 3 xRe{(£ a (s,x.)&(0))} 
2a' G k B T8 al 3, 



(109) 



where we have defined a' G = A 2 /16. The above formula 
shows that a fluctuation-dissipation theorem holds. This 
should be expected because when we expanded the 
determinant and kept only second order terms we were, 
precisely, doing linear response theory. 



III. COMPOSITE MODEL 

We now collect the results from sections U and [TT] in 
a theory modelling the interaction of a spin vector field 
with phonons and electrons. It is straightforward to see 
that the spin-phonon interaction of ^ generalizes to 

m = ~J2(H*^(k)dl(k)S a (k) 



+H a ^k)S a (-k)a fl (k)), 



(110) 



in order to account for a spin vector field (each space po- 
sition now has an independent spin-j degree of freedom) . 
The total Hamiltonian is now given by the Hamiltonian 
of section [il] plus the phonon reservoir Hamiltonian and 
this interaction Hamiltonian. 

Assuming again the factorization of the initial density 
matrix, the Feynman- Vernon functional of this model 
factorizes into the product of two functionals, one as- 
sociated to the phonons and the other to the electrons, 

J^[S l! S 2 ]=J- J) [S 1) S 2 ]J- e [S 1 ,S 2 ]. (Ill) 

The functional associated with the electrons is approxi- 



mated by the exponential of the effective action of (96) 
and the one associated with the phonons is given by 



•F P [S,D] = exp|-i^ dt x J^ dt 2 y / (V2S a (h,-k) ^A«(*i, -k)) G« 7 (ii - * 2 ,k) ^ 



2S 7 (t 2 ,k) 
^D y (t 2 ,k) 



(112) 



where, 



iG ai (t,k)=J2K^)< 



— iw M (k)t 



( -9(-t) 



I 0(f) l + 2n(w„(k)) 



-(A) a7 (i,k)0(-t) 
(A) Q7 (i,k)0(i) 2(A j9 )a 7 (t,k) 



(113) 



r 



Regarding the introduction of the random field in the the sum of two terms. The one arising from the electrons 
equations of motion, now its correlations are given by 
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has been calculated in section [TTJ The contribution given 
by the phonons is given by 



<£„(*, x)$8(0)> 



phonons 



E' 



( i(kx— w M (k)t) 



/r ; ,(k.// ;/! (k, < oti 1 ('^| (k l 



(114) 



Assuming that H afJi (k) = iJ Q (w M (k)), this last expression 
can be written as 



(^(t, X )e /3 (o)) ph( 



dwp(w,x)e-^( k )* 



x J ff*(^(k))iJ )3 (^(k))coth 



in which 



x ) = E S ( w ~ ^W)' 



ik-x 



(115) 



(116) 



If we make the following additional assumption that 

p(w,x) = p(u;)5 3 (x), (117) 

then the discussion of the simplified model of section [I] 
applies to this part of the correlation function and we 
have, for the full correlation function, 

Re{<&(t,x)&(0))} = 6 a0 (a G Lp(t)5 3 (x) 

+ yRe{[J- 1 . 9 ](x)}),(118) 

in which we have applied the considerations we have 
made in section [T] regarding the coupling constants. The 
discussion of the validity of the fluctuation-dissipation 
theorem is the same as in the end of section U and HH It 
is clearly satisfied by both contributions. 



IV. CONCLUSIONS 

Even though the model introduced is very complex, 
we were able to obtain some interesting results. In the 
limit of high temperatures, the magnetic moments associ- 
ated to the 3d electrons feel a random field which has two 
contributions: one from the interaction with the phonons 
which, with some additional assumptions regarding cou- 
pling constants, can be reduced to the one predicted by 
BrowrP^, and the other which comes from the interac- 
tion with the conduction electrons. Besides an exter- 
nal magnetic field, the magnetic moments also feel an 
effective magnetic field associated with their interaction 
with the electrons. This effective magnetic field, given by 
( 131 1, is related to the magnetic field felt by the electrons, 



h(k), and explicitly manifests the Fermi-Dirac statistics 
of these particles because it is weighted by the function 
(Sg)(k) = (l/2)[£ k , (1 - /(e(k'))) /(e(k' - k)) - (k -> 
— k)], where / is the Fermi-Dirac distribution. It is re- 
markable that one can obtain the Landau-Lifshitz-Bloch 



equation from the interaction with phonons using the re- 
sults of Rebei et. aZP and GaranirP. The limit of high 
temperatures should be further investigated in the case of 
the contribution given by the conduction electrons since 
it is not clear from the expression of the random field 
correlation function that it will have Brown's form, i.e. 
white noise, under some assumption of the conduction 
electrons' energy spectrum and density of states in such 
a limit. If this is case, then the Gilbert damping constant 
is given by the sum of two terms, one from the phonons 
and the other from the electrons. The one given by the 
conduction electrons appears to be independent of the 
electron density of states and is equal to A 2 /16. 

The case of finite temperatures is described by two 
generalized functions which are essentially Fourier trans- 
forms of functions which characterize the type of inter- 
action, the statistics and the density of states of the 
bath degrees of freedom. The part of the correlation 
function of the random field which comes from the in- 
teraction with the phonons yields an intimate relation 
between friction (the associated Gilbert constant) and 
the random field fluctuations, see (62). Remarkably, 



this fluctuation-dissipation theorem is manifested even in 
the non-Markovian regime. The validity of the Marko- 
vian approximation is measured in terms of the ther- 
mal correlation time rxh- For time scales lower than 
this, the Markovian approximation for the stochastic 
field fails. In the case of the contribution for the ran- 
dom field correlation function given by the electrons, we 
see that the equation ( 109 ) manifests the existence of a 



fluctuation-dissipation theorem as it should due to the 
expansion done being the associated linear response the- 
ory. The theory considered here satisfies, thus, a general 
fluctuation-dissipation theorem which relates the random 
field fluctuations to the friction constants which measure 
the effect of the interaction of the spin with electrons and 
phonons. 

The model is more general than the so-called three 
temperature model considered and validated in exper- 
imental environment by Beaurepaire et. aZpl This is 
because we do not consider that the spins, electrons and 
phonons are thermalized. Instead, we consider that at 
an initial time to they are thermalized and the density 
matrix decouples at that time, but the time evolution 
couples the various systems and, thus, considering indi- 
vidual temperatures at each instant has no precise mean- 
ing within this model. 

The model proposed here should be further investi- 
gated, namely numerically, since it might give a good 
description of the dynamics of magnetic moments in the 
case of ultrashort laser induced excitations of ferromag- 
netic materials at time scales below the picosecond. 
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APPENDIX: EQUATIONS OF MOTION 

The equations of motion in momentum space of the full 
theory resulting from taking the variation of the effective 



action (in which we consider the truncated expansion of 
|n| appearing in the path integral representation of the 
spin density matrix are 



in which 



'to 



[S(s,k-p) 


x W(s, p) + D(s, k - p) x (T(s, p) 4 


-H eff (s,p))],(119) 


s,p) +H eff 


(s,p)) + *D( S ,k-p)x (jk 2 D(s,p) 


+ W( s ,p))],(120) 


T(s,- 


-k) = T (s) (s,-k) + T (D) (s,-k), 


(121) 




-k) = T^ ) ( S ,-k)+Tif ) ( S ,-k), 


(122) 




-k) = T^( S ,-k)+T^( S ,-k), 


(123) 


' \ r a / 

3 — s) [A[s 


- s , — k) — A{s — s, kjj S(s , — k), 


(124) 


k) = -i\ 2 


f ds'(5 ff )(-k)S(s',-k)6»(s'-s), 


(125) 


■s) [Ap(s- 


s', -k) - A (s' - s, k)] D(s', -k), 


(126) 




= -i^ /^(^(kjD^.-k), 


(127) 


W(s, 


-k) = W ph ( S) -k) + W sd (s,-k), 


(128) 


-s')[A(s- 


- s', -k) - A(s' - s, k)] D(s', -k), 


(129) 


-k) = iA 2 j 


/" rfs'(5g)(-k)D(s',-k)0(s - s'), 

to 


(130) 


= H(-k)H 


-iX [ ds'(Sg)(k)h(-k')6(s' -s), 

J t 


(131) 



r 



here the sub-indices "ph" and "sd" emphasize that this berg exchange term in the action so that it becomes a 

terms come from the interaction with phonon or with diffusive term in the equation of motion in which J is 

electrons, respectively. In the above expression the long some constant measuring the interaction strength, 
wavelength approximation has been taken for the Heisen- 
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